Macroscopic car condensation in a parking garage 
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An asymmetric exclusion process type process, where cars move forward along a closed road that 
starts and terminates at a parking garage, displays dynamic phase transitions into two types of 
condensate phases where the garage becomes macroscopically occupied. The total car density p 
and the exit probability a from the garage are the two control parameters. At the transition, the 
number of parked cars N p diverges in both cases, with the length of the road N s , as N p ~ N% p 
with y p — 1/2. Towards the transition, the number of parked cars vanishes as N p ~ e 13 with /3 = 1, 
e = \a — a* | or e = \p* — p a \, being the distance from the transition. The transition into the normal 
phase represents also the onset of transmission of information through the garage. This gives rise to 
unusual parked car autocorrelations and car density profiles near the garage, which depend strongly 
on the group velocity of the fluctuations along the road. 

PACS number(s): 64.60. Cn, 05.70.Ln, 05.40.-a, 02.50.Ga 
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I. INTRODUCTION 

The scaling properties of nonequilibrium-type phase 
transitions are a topic of intensive current research [1—4]. 
From the theoretical side, focus has been mostly on model 
calculations using stochastic dynamics with simple local 
rules and local interactions. These models serve as pro- 
totypes for various physical processes, such as surface 
catalysis [5], population growth [6], surface growth [2,7], 
electronic transport in wires [8], traffic flow [9], and 
avalanches in granular materials [10]. The simplicity of 
the models is justifiable when, as we expect, the scale 
invariant collective fluctuations at large length scales are 
universal and insensitive to most microscopic details. 

Queuing phenomena, like traffic jams behind slow mov- 
ing trucks on single-lane highways, are common to many 
nonequilibrium processes. One of the interesting phe- 
nomena in such systems is the transition from a finite 
queue to an infinitely long one, i.e., the transition from a 
queue that does not scale with the road length N s to one 
whose length is proportional to N s . Such transitions can 
be induced by changing the total density of cars in the 
system or by varying the probability rate at which cars 
pass the truck [11]. 

A simplification of the above is first to replace the 
slow moving truck by a stationary object, and then to 
ignore the excluded volume effects in the queue by al- 
lowing the cars behind this stationary impurity to pile 
on top of each other. This bare-bone version of the phe- 
nomenon can then be rephrased in terms of the parking 
garage model introduced in this paper. An asymmet- 
ric exclusion process (AEP) on a road starting and ter- 
minating at a garage with an infinite parking capacity. 
The control parameters are the total number of cars, N c , 
in the system (which compared to the road length N s , 
defines a total car density p D — N c /N s ) and the prob- 
ability a with which cars exit the garage (compared to 
the hopping probability along the road, set equal to one). 



The infinitely long queue is represented by a macroscopic 
occupation of the garage. The density of parked cars, 
p,p = Np/N s , is nonzero when the number of cars at the 
garage N p is proportional to N s . 

The analogy with equilibrium Bose condensation is ob- 
vious as well as the basic question, i.e., whether the scal- 
ing properties of such queuing transitions are universal. 
For example, how sensitive are the exponents f3 and x p 
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to the details of the queuing dynamics, particularly to 
the simplifications we made. As a start, we need to es- 
tablish and understand the scaling properties of the most 
simplest model in as much detail as possible. The scal- 
ing properties of the parked car condensation transitions 
presented here are indeed already surprisingly rich. Af- 
ter this, we will be well positioned to restore the omit- 
ted queuing details one-by-one, and establish how robust 
those scaling properties are. 

This paper is organized as follows. The parking garage 
model is introduced in Sec. II, and the phase diagram is 
presented in Sec. III. In Sec. IV, we explore the rela- 
tion to Kardar-Parisi-Zhang (KPZ) type surface growth. 
The group velocity is introduced in Sec. V, and its im- 
portant role in setting the density profile in the various 
phases is reviewed in Sec. VI. In Sec. VII, we do the 
same for the density-density correlations. Sections VIII 
and IX form the core of this paper. We present our nu- 
merical Monte Carlo analysis for the parked car density 
at the condensate-normal phase transitions in Sec. VIII, 
together with explanations of these results. Next, we 
do the same for the parked car fluctuations in various 
phases and at the phase boundaries in Sec. IX. Finally, 
in Sec. X, we summarize our results. 



1 



II. PARKING GARAGE MODEL 

Consider a one-dimensional road of N s sites with peri- 
odic boundary conditions (PBC). Each site on this road, 
1 < x < N s — 1, is either empty or occupied by only one 
car, n x = 0, 1. The parking garage at site N s has no oc- 
cupation upper limit, N p = i%n 3 = 0, 1, 2, 3, .... The total 
number of cars (N c ) in the systems is conserved. From 
the perspective of other physical processes, the cars can 
be reinterpreted just as well as, e.g., molecules (driven 
down a circular tube), kinks in steps on crystal surfaces, 
electrons driven around a wire by an electric field, or 
domain walls in magnetic spin chains. 

The stochastic update rule, from time t to t + 1 is se- 
quential. One of the sites, 1 < x < N s , is selected at 
random with uniform probability p(x) = 1/N S . If that 
site is part of the road, 1 < x < N s — 2, and occupied, the 
car moves forward to over one unit, presuming site x + 1 
is empty; otherwise nothing happens. If the selected site 
is the one just in front of the garage, N s — 1, and is occu- 
pied, then that car moves with certainty into the garage, 
irrespective of the occupation level of the garage. If the 
chosen site is the garage, x — N s , the probability for a car 
to jump out of it onto site x — 1 is equal to a, indepen- 
dent of the number of cars at the garage, provided that 
there is at least one available and site x = 1 is empty. 
This exit probability from the garage is smaller than the 
hopping probability on the road, < a < 1. 

The above process has two control parameters: the to- 
tal density of cars p and the probability a with which 
cars can escape from the garage. We define the following 
quantities: the total car density p — N c /N s , the parked 
car density p p = N p /N s , and the on-the-road density 
Pr = Po — P P - Only p p and p r fluctuate. 

The parking garage is the new aspect to this otherwise 
well-studied AEP. The latter is exactly soluble by the 
Bethe ansatz for periodic boundary conditions (a closed 
loop road without garage) [12], and the stationary state 
properties are exactly known for an open road with reser- 
voirs on both sides [14]. We will use and comment on 
these different setups and exact results in the following 
sections. 



III. THE PHASE DIAGRAM 

Figure 1 shows our phase diagram. It contains three 
phases. In two of them the garage is macroscopically 
occupied: the condensate (C) phase (where the garage 
controls the density of cars on the road), and the max- 
imal current (MC) phase (where the road capacity con- 
trols the flow). In the normal (N) phase, the garage 
contains a finite number of parked cars (typically only 
a few). This structure can be deduced almost com- 
pletely from earlier results for AEPs without the garage, 
such as the AEP with closed periodic boundary condi- 
tions (no garage), i.e., KPZ growth [12,13], and the AEP 



with open boundary conditions hooked-up to reservoirs 
on both ends [14,15]. 

In the two condensate phases, C and MC, the garage 
acts like a reservoir, and the model reduces to the open 
road version studied by Derrida et al. [14] and others [15]. 
In their version, the road is in contact with car reservoirs 
at both ends, x = 1 and x = N s — 1, such that a car jumps 
onto x = 1 with probability a (if empty) and leaves from 
x = N s — 1 with probability (3 (if occupied) . In our model 
(3 is always equal to one. 

A dynamic second-order phase transition takes place 
between the MC phase (where the road density is at a 
constant value, p r = ^) an d the C phase (where the 
road density p r = a varies with the inlet probability a 
and this result is already correct within mean field the- 
ory [14,15]). The MC phase appears because raising the 
density any further would reduce the flow efficiency (due 
to overcrowding). This phenomenon has been canonized 
recently into a "maximal current principle" [16], which 
states that the bulk road density takes the value that 
maximizes the flow. 

Reduce the total number of cars inside the C and MC 
phases, i.e., walk in phase diagram towards the N phase 
along a line of constant a. This has no effect on the 
cars on the road and their fluctuations because all re- 
moved cars are taken from the parked ones residing in- 
side the garage. These surplus cars are dynamically in- 
ert. This continues until the reservoir is depleted, the 
point at which the transition to the N phase takes place. 
The road density becomes equal to the total car density 
p r = p a . From the C (input-limited) phase perspective, 
this happens at p = a, see Fig. 1, because p r = a inside 
the C phase [14-16]. From the MC (road- limited) phase 
perspective, it occurs at p = ^ because p r = ^ inside 
MC. Therefore, the phase boundaries into the N phase 
are located at p Q = a and p a = 5 - 

The next issue is how many types of normal phases 
exist in the lower right side of the phase diagram. In 
a N phase, the garage acts like a localized impurity (a 
blocking- type site), and it contains typically only a few 
cars. Suppose we put a cap N m on the occupation of the 
garage, N p = 0,1, N m . This cannot affect the prop- 
erties of the model in the N phase, except very close 
to the transition into the condensate phases (where the 
number of parked cars N p and the fluctuations in N p di- 
verge). Janowsky and Lebowitz [11] (JL) studied such 
a capped version of the AEP model, a closed loop road 
with periodic boundary conditions, without a garage, but 
with one special bond where the hopping probability is 
reduced from 1 to a. This is equivalent to setting the 
occupation limit of the garage to N m = 1. Their phase 
diagram has the same control parameters as ours. It 
has particle-hole symmetry with respect to p a = |, see 
Fig. 2. The normal phase in the lower right hand cor- 
ner is similar to our model. The second normal phase in 
the upper right corner is equivalent to it by particle-hole 
symmetry, and the intermediate area is a coexistence re- 
gion, the jammed phase. In the low-density N phase, the 
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cars are uniformly distributed, but at p Q = a/(l + a) a 
traffic-jam- type shock wave develops between low- and 
high-density TV- type phase regions. 

This jammed phase is reminiscent of the condensate 
phase (macroscopic occupation of the parking garage) in 
our model. A shock wave does not appear in our model 
because the garage absorbs the queue. It is also clear 
that if we would interpolate between the JL model and 
our model by slowly increasing N m = 1 — ► oo, the queu- 
ing transition would be delayed and would be shifting 
smoothly, with our N-C phase boundary as the limiting 
case. JL focused on the finite-size-scaling (FSS) prop- 
erties of the fluctuations in the position of the tail of 
the shock wave, Al <~ L 1 , and found 7 = | (except for 
p = 2 where a cancellation of the leading time of flight- 
type fluctuations gives rise to 7 = 3) [11]. Monitoring the 
location of the shock wave is analogous to observing the 
parked car density fluctuations in our model inside the C 
and MC phases. JL did not address the scaling proper- 
ties at the queuing transition itself, but their transition 
is definitely second order (the location of the shock wave 
moves continuously with the total density away from the 
blockage point). 



IV. KPZ GROWTH 

The AEP can be interpreted as a lattice representation 
of the Burgers equation with a discretized velocity field, 
both in location and in its values, n(x) = 0, 1. The latter 
is equivalent to the so-called body-centered solid-on-solid 
growth model [7,17], a lattice realization of KPZ-type 
surface growth. Each occupied site represents a down 
step along the one-dimentional (ID) surface and every 
empty site an up step. The surface heights, h(x + i), 
are limited to even/odd values at even/odd lattice sites, 
h(x+7j) — h(x— \) = 1 — 2n(x). The surface looks like the 
alternatively stacked bricks of a masonary wall, but ro- 
tated over 90 deg. Each car jump to the right represents 
the deposition of a new vertical 1x2 brick. The average 
car density on the road represents the average slope of 
the surface. At half density the surface is not tilted. 

Periodic boundary conditions are the most common in 
theoretical studies of surface growth, and this brick de- 
position version of KPZ growth is exactly soluble by the 
Bethe ansatz [12]. It is a special line in the general phase 
diagram of the ID well-known XXZ spin-^ quantum spin 
chain Hamiltonian, and also in the so-called six-vertex 
model of 2D equilibrium critical phenomena [18]. 

From the surface growth perspective, the open road 
setup with reservoirs represents an open ended ID in- 
terface with modified particle deposition probabilities a 
and (3 at its end points. In the MC phase, where the car 
density is locked to p a = | , the crystal surface maintains 
a net zero tilt. For a < i and/or (3 < |, the reduced 
growth probabilities at the edge create a nonzero globally 
tilted surface. Moreover, along the line a — [3 < \, two 



tilts coexist. Crossing that line amounts to undergoing a 
first-order phase transition with different tilt angles. The 
lock-in transition from the C phase into the MC phase 
represents a second-order dynamic facetting transition. 

Our parking garage setup translates into a surface 
growth dynamics realization with periodic boundary con- 
ditions and a localized defect. The garage could represent 
something like a stacking defect, where the step height is 
unlimited (and somehow energetically cost free). The cliff 
height at the stacking fault can be microscopic (the N 
phase) or macroscopic (the C and MC phases). This de- 
pends on the modified growth probability a near the de- 
fect (on top of the cliff, somewhat similar to a Schwoebel 
energy barrier) and also depends on the global net tilt 
angle (the total car density). In the MC (C) phase, the 
average slope of the surface, excluding the cliff, is zero 
(nonzero). 



V. CORRELATIONS AND GROUP VELOCITY 

It is well known and easy to show that the station- 
ary state of KPZ growth with PBC, the closed loop road 
without obstacles where a = 1 and n(N s ) = 0, 1, is com- 
pletely random without any correlations whatsoever be- 
tween finding a car at sites i and j for any i ^ j including 
nearest neighbor sites. In the surface representation, this 
means that up and down steps are placed at random [17] 
and that the width of the interface W ~ N* scales with 
the random walk exponent X—\- 

The stationary state remains uncorrelated in our model 
as well (as shown below) except near the garage. Those 
correlations are governed by the dynamic scaling ex- 
ponent and the group velocity with which fluctuations 
travel along the road. 

In ID KPZ-type surface growth, all characteristic 
times scale as t <~ N* with z = 3/2. In car language, 
this means that local car density fluctuations with spa- 
tial width I broaden in time as I ~ t 1 ^. The value 
of the dynamic exponent z follows from the identity 
z + x = 2 implied by Galilean invariance of the Burg- 
ers equation [12,13], together with the disordered x — \ 
nature of the stationary state. The exact (Bethe ansatz) 
solution of this specific model confirms this general re- 
sult. 

Fluctuations travel with a group velocity v g = 1 — 2p r 
to the right (towards the garage) along the road. From 
the KPZ surface growth perspective, the group velocity 
represents the average slope of the KPZ growing surface, 
and its fluctuations grow perpendicular to the surface 
orientation. From the AEP formulation perspective, the 
value of v g is set by the average stationary state current 
j = p r (l — p r ) (easy to derive since the stationary state 
is Gaussian) and the definition of the group velocity is 

d 

V 9 = Ttrj(Pr) = (! - 2 Pr)- ( 3 ) 



3 



This v g played an important role, e.g., in the AEP study 
by Majumdar and Barma [19] of tagged-particle diffu- 
sion, and in describing phase transitions between steady 
states in the open road setup [20] . 



VI. DENSITY PROFILES IN THE VARIOUS 
PHASES 

The on-the-road car density profiles in the C and MC 
phases are known exactly from the open boundary model 
studies, in particular the one by Derrida et al. [14]. These 
profiles play an important role in our later discussions 
and it is useful to review how the group velocity enters 
in qualitative explanations of these exact results. 

In the road- limited MC phase at a > \, the density 
tail has a power-law shape, 



p(x) ~ p r + A e X 



(4) 



at both edges e=L(left),i?(right), with x the distance 
from the edge of the road (garage entrance in our model) 
and p r = \ the on-the-road car density. In contrast, 
in the input-limited C phase at a < \, the car den- 
sity is constant, p(x) = p r , at the beginning of the road 
(garage exit), and has an exponential tail at the end of 
road (garage entrance), 



p(N s — x) ~ p r + Ax 



-3/2 e -x/« 



(5) 



with the on-the-road car density p r — a and the correla- 
tion length £ = -l/ln[4a(l - a)]. 

Power laws with exponent v — ^ arise naturally in this 
problem because of its critical fluctuations. The bulk 
properties of the KPZ stationary state are invariant un- 
der a rescaling of all lengths as x' — bx, all times as 
t' = b z t, and the surface heights as h' — b x h. The car 
density scales, therefore, naively as p — dh/dx ~ 6 X_1 . 
This means that power-law tails in the density distribu- 
tion near the edge of the road, with exponent \~ 1 = — \, 
like in Eq. (4), are to be expected. 

On the other hand, the disordered nature of the sta- 
tionary state suggests exponential profiles, like in Eq. (5). 
Indeed, in the bulk, the stationary state density-density 
correlation function 



g(r) = (p(x + r)p(x)) - (p(x + r)> (p(x)) 



(6) 



does not decay as a r _1 power law as may be suggested 
by the above scaling argument, but instead decays ex- 
ponentially with a very short correlation length, that is 
zero in this specific model since the cars on the road are 
totally uncorrelated in the stationary state. 

To make sense of Eqs. (4) and (5), it is important to 
realize that the density profile near the edge of the road 
incorporates temporal correlations, and also to appreci- 
ate the role of the group velocity. 

The power-law profile in the MC phase is the result of 
correlations with cars that reached the edge of the road 



and moved into the garage at earlier times. Such corre- 
lations spread in time over a spatial distance I ~ t 1 / 2 . 
The group velocity is zero in the MC phase and there- 
fore the cars at a distance I from the edge of the road are 
(power-law) correlated with those that entered the garage 
at times earlier than At = P. They have no knowledge 
about cars arriving at the garage more recently. This 
explains the power-law tail in the density distribution. 

In the C phase, the same correlation spreading takes 
place within a moving frame of reference with nonzero 
group velocity. Information reaches the edge of the road 
(garage entrance) at a rate v g (linear in time) faster than 
it can spread backwards as I ~ t x l z (since 1/z = 2/3 < 
1). Memory has no opportunity to develop near the edge 
of the road, and the density profile adjusts itself at the 
road end as if the cars on the road are completely un- 
correlated. The exponential tail reflects a suction-type 
effect; the excluded volume limitation on the car mobil- 
ity does not act on the cars departing from the road. 

Similarly, in the MC phase, power-law correlations 
with earlier cars that escaped from the garage and en- 
tered the road give rise to a power-law tail in the density 
profile at the beginning of the road (near the exit of the 
garage), while in the C phase, this information travels 
faster away from the road start than it can spread back- 
wards. So, new cars entering the road do so completely 
uncorrelated, resulting in no road-start tail in the profile 
whatsoever. This confirms why the MC phase is road 
limited and the C phase is input limited. In the latter, 
the car supply at the garage controls the density near the 
edge of the road and also everywhere else on the road. 

In the TV phase, we find numerically, from Monte Carlo 
simulations, an exponential tail at the road end and an 
1/x tail at the road start as shown in Fig. 3. The group 
velocity is also nonzero in the iV phase, which explains 
the exponential exit tail as follows. Just like in the C 
phase, the garage is invisible to the incoming of cars 
from the road; fluctuations arrive at the garage faster 
(at constant velocity) than they can spread backwards 
(with Al <~ i 1 / 2 ). The fact that in the TV phase only a 
few cars reside inside the garage is invisible to the cars 
entering the garage. 

An 1/x tail at the beginning of the road is what one 
expects from the deterministic part of the Burgers equa- 
tion. The solution of the deterministic Burgers equation, 



dv dv d 2 v 
dt dx dx 2 



(7) 



with the velocity v pinned at a specific nonzero value at 
site x — 0, is readily seen (e.g., in the Hopf transformed 
formulation) to be of the generic form v(x) <~ 1/(1 + ax), 
i.e., having the 1/x power- law shape. 

In the C phase, the system selects the value a = for 
the constant, and in the N phase, s/0. In both cases, 
the group velocity v g > carries away KPZ fluctuations 
faster than they can spread, such that the KPZ noise 
term can be ignored at the road start (stationary frame 
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of reference) . The only noise that remains is that from 
the random process by which cars are being taken out 
of the garage and put on the road. In the C phase, the 
supply of cars is bottomless (p p > 0), such that a = 
(the entire on the on-the-road bulk car density is ruled 
by the garage), while, in the N phase, the supply of cars 
is limited and the garage does not overwhelm the road, 
such that a ^ 0. 

VII. DENSITY-DENSITY CORRELATIONS 

It is useful to discuss how the group velocity affects 
the car-car correlation function, 

g(r, t) = (p(x + r) t+T p(x) t ) - (p(x + r) t+T ) (p(x) t ) . (8) 

According to scaling theory, it obeys the form 

g (r,T) = b 2 ^- 1 ^g(b- 1 r,b- z T) 
2(v— i) r 
= -^n— z )- 0) 

In the limit (f> = t/t 1 / 2 — > 0, the scaling function F{4>) 
must approach a constant since the autocorrelation de- 
cays as a power law, 

. 9 (0,t)~t^~t- 2 / 3 , (10) 

with x = \- ln the opposite limit, of large (j>, the scaling 
function must decay exponentially because g(r, 0) = 
due to the random nature of the stationary state (or, 
more generically, it decays exponentially with some cor- 
relation range 4> of the same order as the interaction 
range between cars) . The above scaling relation suggests 
the form 

F(4>) ~ 4> 2 ^) e -<t>/<t>a (11) 

in the limit <j> = r/r 1 / 2 — > oo, such that 

g(r, t) ~ r 2 (x-i) e -^^ 1A ^ l e -r/i ( 12 ) 

with x = 1/2 and I = (^qT 1 ^. The length I in the ex- 
ponential defines the "correlation cone" that we already 
mentioned above. Within the cone the correlations decay 
as a power law, and outside it the cars are uncorrelated 
(as in the stationary state). The above discussion ignores 
finite-size effects. For open road or garage type bound- 
ary conditions, the correlator explicitly depends on the 
initial position xq of the car at time to and also on the 
road length N s , 

g(r, r; x , N s ) = b 2 ^ g{b~\ b~ z r; b^xo^N.). 

(13) 

Finite-size effects set in at times when the correlation 
cones hit the road edges, r ~ (xo) z or t ~ (N s — x n ) z . 



At times longer than to <~ N z , we therefore expect that 
the car-car correlator decays exponentially, e.g., as 

g (o,t; lN s ,N s ^j ~ N 2 ^e'^ N ^. (14) 

The above analysis applies to the MC phase, where the 
drift (group) velocity is zero. In the C phase (and also the 
N phase), we need to switch to the moving frame of ref- 
erence by replacing r — > r + v g r. This has some peculiar 
consequences. For example, the autocorrelation function 
g(0, t) decays exponentially in time (even at short times). 
The correlation cone I is slanted in the direction of the 
flow (the correlations move with the flow) and I widens 
slower than linear in time, such that the r = line in 
the world sheet lies outside the correlation cone. The 
KPZ-typc correlations are somewhat hidden. In order 
to expose them, one needs to plot them in some special 
manner, e.g., g(v g T,T) (the autocorrelator in the moving 
frame of reference), as illustrated in Fig. 4. 

VIII. SCALING AT AND NEAR THE 
CONDENSATION TRANSITIONS 

Let us turn now to the scaling properties of the two 
condensate phase transitions from the C and MC phases 
into the N phase. The stationary state value of the 
parked car density p p — N p /N s acts as the order pa- 
rameter. We expect it to obey the FSS relation 

p p (e,N- 1 ) = b-^p p (by^e,bN- 1 ) (15) 

near the two condensation transition lines with a scaling 
factor b, and e = (a — a*) or e = (p* — p a ), a measure of 
the distance from the transition. 

In the TV phase, the density of parked cars is zero, 
p p = 0. When approached from the C and MC sides, 
p p goes to zero as p p <~ |e|*, where (3 = x p /y t . The 
removal of cars from the road is accommodated by tak- 
ing them from the passive inert ones residing inside the 
garage. This implies that p p must vanishes linearly at the 
transition, and that x p — y e . Our numerical simulations 
confirm that p p <~ |e|, as [3 = 1 in Figs. 5 and 6. 

The total number of parked cars (N p ) scales as 

NpicN- 1 ) = by-Npiby-e^N- 1 ) (16) 

with exponent y p — 1 — x p according to Eq. (16) and 
N p = p p N s . In the TV phase, the density p p remains 
zero, while the total number of parked cars (N p ) di- 
verges towards the transition as N p ~ \e\ < - 1 ^ x p^ Vc . We 
find numerically that this power law is linear as well, 
N p ~ |e|. Combined with the linear scaling of the car 
density from the C or MC side (implying y p = x p ), this 
yields y p — y t — \. The exponent y p determines the FSS 
behavior of the total number of parked cars at the tran- 
sition point itself, N p <~ iVf p (and also that the density 
of parked cars vanishes as p p <~ 
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Figure 5 shows the numerical results of y p and (3 at 
point a — 0.25 of the C-N phase boundary and Fig. 6 at 
a = 1 (and 0.75) of the MC-iV phase boundary. Notice 
that the FSS corrections in the exponent y p are much 
stronger at the MC-iV transition than at the C-N tran- 
sition. 

The validity of the scaling relations is further illus- 
trated by plotting the scaling function <!>(£) in the fol- 
lowing form 

p p (e,N- 1 ) = N- x ^(N^e), (17) 

and the associated data collapse while moving through 
the transition points at constant p Q and constant a, re- 
spectively. The curves in Fig. 7 collapse very well. 

The fluctuations and the FSS corrections to the num- 
ber of parked cars are a mirror of the density-density 
correlations and the density profile of cars on the road. 
They also reflect how the latter builds up as a function of 
the length of the road. At the transition points, the bulk 
value of the on-the-road density p r = N r /N s is exactly 
equal to the total number of cars in the system divided 
by the road length, p r = N c /N s , such that there would 
be no need for any car to remain inside the garage. 

We find numerically that at the C-N transition point, 
the density profile retains the same structure as inside 
the C phase; with no tail at the beginning of the road 
(garage exit) and an exponential tail at the end of the 
road (garage entrance) . Such profiles cannot account for 
the y p = \ FSS divergence in the number of parked cars. 
The scaling behavior N p <~ Ns p must, therefore, reflect 
directly the corrections to FSS in the bulk density of cars 
on the road. The value x p — | naturally arises because 
at the transition point, p p ~ p r , and p r = dh/dx scales 
as L x_1 ~ L~ Xp , where L corresponds to a given length 
of the road N s . 

A more intuitive explanation follows again from the 
nonzero group velocity at the C-N phase boundary. As 
mentioned above in the discussion of the density distri- 
bution tail in the C phase, the events by which cars enter 
the road from the garage are completely uncorrelated due 
to the slanting of the correlation cones (the correlations 
move with the flow towards the garage and spread slower 
than linearly, only as I ~ t 1 ^, such that communications 
with later events at the beginning of the road are impos- 
sible). So the entry events to the road from the garage 
behave like uncorrelated random noise, and the fluctua- 
tions in the number of cars scale, therefore, as the square 
root of time. The time in question is proportional to 
N s because the fluctuations (created at the garage exit) 
move along the road with velocity v g , and are wiped out 
after they return to the garage. From this, it follows that 
the fluctuations in the number of parked cars scales as 
VN S . Moreover, N p cannot be negative, which means 
that the fluctuations sample the bottom of the garage 
and that the transition from the C to N phase, there- 
fore, takes place when the garage contains N p <~ VN S 
cars. 



In summary, the scaling at the C-N phase boundary 
is governed by the bulk fluctuations in the-on-the-road 
density, which is ruled by the nonzero group velocity of 
KPZ fluctuations, and this leads directly to random noise 
like N p ~ \/N a , corrections to FSS. 

The FSS corrections to scaling in exponent x p (y p ) at 
the MC-7V phase boundary are much more complex. The 
(bulk) group velocity is zero, and power-law density pro- 
files are realized at both edges of the road. At the end of 
the road, the density profile follows a critical exponent 
v = 5, the same power as that inside the MC phase dis- 
cussed in Sec. VI. However, at the road start (the exit of 
the garage), the power-law exponent changes from v = \ 
inside the MC phase to v = | at the MC-7V transition. 
This is shown in Fig. 8 for a = 1 and p a = ^, where the 
effect is the strongest. 

The I power law does not change the Nl^ 2 FSS be- 
havior of the number of parked cars. It is responsible, 
however, for strong corrections to FSS, as clearly visible 
in Fig. 6(a). The v = | power-law profile contributes 
only a subdominant term to p p because it decays faster 
than the two v = 5 contributions (from the density pro- 
file at the end of the road and from the KPZ-like bulk 
road density fluctuations). We have not achieved yet a 
good understanding of this novel value, |, for the expo- 
nent of the density profile at the beginning of the road. 
It obviously lies correctly in between the MC and N val- 
ues, ^ and 1, respectively. Moreover, its value, v = |, 
is likely linked to the KPZ dynamic exponent z = |. 
But it remains unclear how to glue the following pieces 
together. 

At the MC-iV transition, the group velocity of the fluc- 
tuations, v g = 1 — 2p r , is still zero (but only barely) 
since p r = |. This means that time of flight aspects, 
which dominate the C and N phases, do not come into 
play. Near the garage, the tail in the density profile, 
p{x) = I + Ap(x), creates a local nonzero group velocity 
v g (x) = —2Ap(x) pointing back into the garage (instead 
of away from it as in the C and N phases). 

Recall from Sec. VI that the density profile in the MC 
phase has the same type of power-law profile Ap(x) ~ 
x~", but with the "KPZ" (power-counting) exponent 
value v = |. In that case, fluctuations at the road start 
do not reach the bulk of the road: the leading edge of 
the information cone is stationary because the backward 
movement of its center of mass x c ~ — (implied 
by dxcjdt <~ x~ v ) matches exactly the rate of its spread- 
ing, x ~ i 1 / 2 . So the v = \ density profile fully screens 
the garage from view in the bulk of the road. Exactly 
the same screening of the information takes place at the 
end of the road (garage entrance) since also there the 
MC phase density profile has exponent v = 5 (but in an 
opposite forward moving v g sense and with a negative 
density profile amplitude). 

The density profile exponent v = | at the MC-iV tran- 
sition does not fully screen the garage any more from 
observers located far away on the road. Total screening 



G 



is not needed because KPZ fluctuations start to tunnel 
through the garage since it has only a y/N s occupation. 
It is yet unclear to us, however, how to deduce (in a 
convincing manner) the exponent v = | from these con- 
siderations. 



IX. PARKED CAR FLUCTUATIONS 

In this section, we explore the fluctuations in the 
parked car density and also car-car correlations on the 
road. Of particular interest is the onset of transmission 
of information through the garage at the phase transi- 
tions. In the two condensate phases, the garage acts as 
a car reservoir and as a sink of fluctuations, while in the 
normal phase it contains only a few cars and transmits 
fluctuations. 

The temporal fluctuations in the total number of 
parked cars, G(N s ,t), measures also the fluctuations in 
the total number of cars on the road. It is therefore equal 
to the integrated car-car correlator (defined in Sec. VII), 

N s -1 

G(N s ,t)= [{p( x o,t )p(xi,h)) -p(x Q )p(xi)] 

Xl,XQ — l 

N s -1 N s -x -l 

= E E 9(r,r;x ,N s ) (18) 

x = l r= — xq + 1 

with r = ti — t a and r = X\ — x . The summations run 
over all road sites xo and distances r that fit on the road. 
Direct integration of the scaling relation, Eq. (13), yields 
that (in the MC phase) G obeys the scaling form 

G(N S , t) = b 2 * G{b- X N S , b~ z r) = N 2x T{t/N z ). (19) 

In the KPZ representation, G is the global slope-slope 
autocorrelator, and at r = 0, reduces to the conven- 
tional definition of interface width (second moment of 
the height distribution). We will now discuss how G be- 
haves in the various phases and at the phase transition 
points. 

A. Inside the MC phase 

The following intuitive discussion tells us how the scal- 
ing function T behaves in the MC phase. Each g(r 7 r) 
has a correlation cone of size I ~ t x I z . The cars within 
this cone are correlated with the one at site xq at time to 
as t 2 (* -1 )/ 2 . The integration over xo and r in Eq. (18) 
yields 

G ~ N s x I x r^x-i)/* ~ A^ 2 *" 1 ^, (20) 

where the correlation cones I <~ r 1 / 2 are assumed to be 
small with respect to the road size N s . The exponent 
is equal to \ = 5 , such that G ~ N s , and that T ap- 
proaches a constant in the limit t/N z — > 0. 



This estimate fails to take into account finite-size ef- 
fects. Correlations are truncated near the two road edges 
(all information is entering the garage). The loss term is 
of the order 

2 [ dx (I - xo)t 2 ^-^' z ~ I 2 l 2{ *- l) ~ r 2 */ z . (21) 
Jo 

This suggests that G is of the form 

G = N.(a-bj£- + --?), (22) 

with constants a and b, and suggests that the scaling 
function F(4>) is analytic at short times in the parame- 
ter <j) = t x I z /N s instead of <j>' = t/N z . Our numerical 
results shown in Fig. 9(a) are consistent with this. 

In the opposite limit, (j>, <fi' — ► oo, where time is large 
compared to the length of the road, all correlation cones 
are limited and equal to I ~ N s , and G behaves as in 
Eq. (14), such that 

G ~ N s x I x g(0, t; ^N s ,N s ) ~ N s e -< T l N '\ (23) 

Our numerical results in Fig. 9(b) are consistent with this 
as well. 



B. Inside the C phase 



In the C phase, the fluctuations scale with the same 
exponents as in the MC phase, but the nonzero group ve- 
locity v g = l — 2p r changes completely the appearance of 
correlation functions, like G(N S , r) . The density-density 
correlations spread just like in the MC phase, but only 
with respect to the moving frame of reference and with r 
replaced by f = r + v g T. The correlation function g(r, t) 
scales as a power law g <~ t 2 ^ 1 ^ 2 at f = v g r (i.e., 
r = 0) but exponentially at nonzero r. The fluctuation 
cones, I ~ i 1 / 2 , are slanted and move with the flow. 

Figure 10 shows G(N s ,t) inside the C phase as a func- 
tion of time t for various road sizes N s . It decays linearly 
until hitting zero at Tfl; g ht = N s /v g (= T), and then it re- 
mains at zero. All fluctuations move with the group ve- 
locity v g to the right and reach the garage at a constant 
rate. After one time of flight, T, all correlations with the 
initial configuration have disappeared. The rounding in 
G at T is of order At ~ T 1 / 2 , and is due to the broad- 
ening of the remaining correlation cones just before they 
are absorbed by the garage. 
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C. Inside the N phase 

In the TV phase, the fluctuations in the total number of 
cars on the road, G(N S , 0), are not proportional to N s but 
are only of order one. The garage is not macroscopically 
occupied any more and acts very much like an ordinary 
road site. Figure 11 shows the behavior of G(N s ,t) in- 
side the N phase. The total number of cars in the system 
is conserved, such that G reduces to 

G(N S , t) = (N p (t + r)N p (t Q )) - N 2 p , (24) 

and behaves similar to the autocorrelator g(0,r). The 
group velocity is nonzero, and therefore G decays expo- 
nentially fast. However, because of PBC, G comes back 
to live, like a lighthouse light beam, after every time-of- 
flight interval T(= N s /v g ), with an amplitude of order 
T 2(x-i)/z and with a tempora i width f orc jer I ~ T 1 ^. 



D. At the C-N transition 

Figure 12 shows how G(N s ,t) at the C-N transition 
decays in time for various system sizes. At small t, it de- 
cays linearly, rather like in the C phase, but then it seems 
to oscillate with a period determined by the time-of-flight 
time scale T; G(N S , r) goes through zero at about t ~ 
and shows a strong anticorrelation at t ~ T (the maxi- 
mum lies just before it). 

Figure 13 illustrates how the transmission of informa- 
tion through the garage commences at the phase tran- 
sition. Suppose we approach the transition point from 
the C phase following a line of constant a. In the C 
phase, these lines coincide with lines of constant group 
velocity. (In the N phase, v g is constant along lines of 
constant p a .) So nothing changes on the road until we 
hit the transition point, and G decays linearly to zero 
and remains zero after one time-of-flight time scale. At 
the transition point, G transforms abruptly into the os- 
cillatory shape, with an anticorrelation after one time of 
flight. After that, it reduces inside the N phase to the 
lighthouse shape, in which G(N s ,0) oscillates in phase 
and does not scale with N s any more. 

The anticorrelations at the C-N transition point are 
intriguing and need to be explained. Imagine a localized 
positive density fluctuation at the beginning of the road 
at time t - In the MC phase, it simply sits there while 
broadening as I ~ t x I z ~ t 2 / 3 and weakening in ampli- 
tude as t 2 ^ 1 ^ 7 - ~ r~ 2 / 3 . In the C phase, it broadens 
and weakens in the same manner, but travels like a soli- 
tary wave to the right with velocity v g = 1 — 2p r and 
drops out of the road after one time-of-flight unit T. In 
the TV phase, it behaves very much the same, except that 
the positive density fluctuation creates a deficit inside 
the garage since the number of cars in the garage is fi- 
nite, such that fewer cars can be put on the road in the 
immediate wake of the positive fluctuation. Therefore, in 



the TV phase, every positive fluctuation carries a compen- 
sating (again localized) negative tail with it. Figure 14 
illustrates the difference schematically, and our numerical 
simulations confirm this picture. 

One could say that in the N phase positive and neg- 
ative local excitations are bound in pairs, and that they 
unbind at the N-C transition. On approach of the tran- 
sition from the N side, the width of the negative tail 
grows, but with conserved total area (equal to the area 
of the positive part of the excitation), because the aver- 
age number of cars in the garage increases towards the 
transition (and diverges) and therefore the reduced car 
output is being spread over more time. At the transition 
point itself, the negative tail has vanished, except for the 
finite-size scaling effect of order (A s ) -1 / 2 . 

Local excitations, thus, behave quite interestingly. 
However, they do not explain the difference in the behav- 
ior of G at the transition point and in the C phase. These 
solitary waves do not deplete the garage sufficiently to 
trigger its bottom because at the transition point the 
number of parked cars still diverges as y/N s . 

Only nonlocal excitations, which encompass the entire 
system, are able to empty out the garage. Consider an 
excitation where the density of cars is globally and uni- 
formly enhanced along the entire road, pi = p r + A. This 
requires that the number of cars to be taken out of the 
garage should be proportional to the road length N s . In 
the C phase, those will not deplete the garage because 
the number of parked cars is also proportional to N s . 
Such a A regiment of extra cars marches with group ve- 
locity v g to the right, reaches the end of the road, and 
thus returns to the garage at a rate uniformly in time, 
row by row. Throughout this process, the garage is not 
aware of the existence of the regiment since we did not hit 
its bottom and because the information cones I ~ r 1 / 2 
on the road do not broaden fast enough compared to the 
group velocity to maintain its memory at the garage exit. 
This implies that the exact stationary state rebuilds itself 
in the wake of the regiment, and also that the enhanced 
road density decays linearly in time and vanishes com- 
pletely after one time of flight T, just like our G as a 
function of time r. 

At the transition point, the garage only contains N p ~ 
VN S cars, so that the same type of global uniform exci- 
tation can only have an amplitude of order (Ns)' 1 / 2 , and 
always depletes the garage. This regiment of cars travels 
to the right with velocity v g , just like in the C phase, 
but in its wake the garage cannot rebuild the station- 
ary state because it is empty. A depleted road density is 
established in the wake of the enhanced excitation, and 
thus anticorrelations build up, and after one time of flight 
the average density of cars on the road is below normal 
(and a surplus of cars resides in the garage). This ex- 
plains qualitatively the oscillatory behavior of G at the 
transition point and the anticorrelations at T. 
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E. At the MC-N transition 

The correlations at the MC-N transition are less spec- 
tacular than at the C-N transition. The group velocity 
is zero inside the MC phase, and still remains zero at the 
MC-N transition. Only inside the N phase, does it start 
to shift continuously away from zero. Figure 15 shows 
how G scales at the MC — N transition. These numer- 
ical results are almost the same as those inside the MC 
phase in Fig. 9. 

Recall the argument about the behavior of G(t, N s ) in- 
side the MC phase (in section IX A), and imagine how 
this was modified at the MC-N transition. The factor 
N s in Eq. (20) represents the number of cars on the road 
(the number of sites x that are occupied). This should 
be modified to p N s — a^f~N s , since the number of parked 
cars scales as VN S . The other terms, the spreading in 
time of the correlation cones and the autocorrelations on 
the road, are likely unchanged. Such differences are sub- 
tle and not surprisingly numerically invisible. 



There is some evidence suggesting that the expo- 
nents are indeed robust, e.g., in the Janowsky and 
Lebowitz [11] model, the same simple KPZ-type values 
of the exponents appear in the fluctuations of the queue, 
although a detailed study of the queuing transition it- 
self needs still to be performed. In addition, introducing 
short-range car-car interactions [15,21] do not seem to 
change the exponents either [22]. 
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X. SUMMARY AND CONCLUSION 

In this paper, we presented the scaling properties of 
dynamic condensate phase transitions in terms of an 
ID asymmetric exclusion process with a parking garage. 
There are two types of condensate phases: the maximal 
current (MC) phase, where the road controls the den- 
sity of cars on the road, and the condensate (C) phase, 
where the garage (as a reservoir) controls the number of 
cars on the road. The existence of a group velocity is cru- 
cial for understanding the behaviors of correlations and 
the density profiles in these two phases and at the phase 
transitions. 

At both condensate transitions, the number of parked 
cars scales as N p <~ Ns" with y p = ^, while on ap- 
proach of the transition, the density of parked cars van- 
ishes linearly with the control parameters (the total den- 
sity of cars in the system and the exit probability from 
the garage) p p <~ \e\@ , with (3=1. Also, the transition 
points represent the onset of communication of informa- 
tion through the garage. This leads to interesting auto- 
correlations in the number of parked cars, particularly 
at the C-N transition, due to the nonzero group velocity 
and associated time-of-flight effects. 

Our parking garage model is a bare-bone version of 
dynamic Bose condensation and of queuing phenomena, 
like traffic jams. The fundamental issue that needs fur- 
ther study is whether the above scaling behavior, in par- 
ticular the values of the critical exponents, are universal 
or not. How do more realistic interactions between the 
cars on the road change this? Do the simplifications in 
the traffic jams, like a stationary truck versus a moving 
one and ignoring the spatial structure inside the queue of 
cars behind it by collapsing them (piling them up) into 
a "garage," affect the exponents? 
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FIG. 2. Phase diagram of the Janowsky-Lebowitz model, 
an AEP with periodic boundary conditions and a slow bond. 
The corresponding density profiles are schematically shown in 
the insets. 




x 

FIG. 3. Double logarithmic plots of the density profiles in 
the N phase (p — 0.25 and a = 1), showing al/i tail at the 
beginning of the road (garage exit). 
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FIG. 1. Phase diagram with the total car density p and 
the exit probability a from the garage. At the phase transi- 
tion lines into the N phase, the parked car density p p vanishes, 
and the garage seizes to be macroscopically occupied (shown 
schematically in the insets). 
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FIG. 4. Double logarithmic plots of the density-density cor- 
relation function g(r,r) in the moving frame of reference, i.e., 



with 



„t and v„ 



(1 



2p r ). The C phase data are 
obtained at p — 0.75 and a = 0.25, the C-N point data at 
p o =a=0.25, and the N phase data at p — 0.25 and a = 1. 



10 




(b)N s =512 
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FIG. 5. The parked car (order parameter) exponents 
(a) y p = 0.505(5) and (b) [3 = 1.000(1), denned as N p ~ N%" 
and p p ~ (a — a*f , at p o —a=0.25 of the C-N phase bound- 
ary. 
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FIG. 6. The parked car (order parameter) exponents 
(a) y p = 0.495(5) and (b) (i = 0.997(3), at p a = 1/2 and 
a = 1 of the N-MC phase boundary. For other two curves 
in (a), taken along the MC-N phase boundary at a > 1/2, 
the exponent y p retains the same value, but is subject to 
large/small corrections to finite-size scaling originating from 
2/3 power- law density profile at the beginning of the road 
(garage exit). 
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FIG. 7. Scaling function defined in Eq. (17) at 

the same transition points as in the previous two figures; 
(a) p = a = 0.25 (C-N) and (b) Po = 1/2, a = 1 (MC-N). 
The numerical data collapse very well for x p = y,_ — 1/2, 
as suggested in Eq. (17). The scaling function $(£) — ► for 
£ « and $(£) -> £ for C > 0. 
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FIG. 8. Effective (finite-size-scaling type) value for the ex- 
ponent v of the density profiles at the beginning of the road: 
v = 0.66(1) at p = 1/2 and a = 1 (or 0.75). At the MC-iV 
phase boundary, the density profiles decay as a 2/3 power law, 
instead of the 1/2 power found inside the MC phase. 
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FIG. 9. The parked car correlator G(N s ,t) in the MC 
phase, for small r (a), it decays linearly as function of 
r 1,/z /N s , and for large r (b), it exponentially as function of 
t/N*, where z = 3/2. The data are obtained at p = 0.75 
and a = 1. 
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FIG. 10. The parked car correlator G(N s ,t) in the C 
phase; it decays linearly as function of r and becomes zero 
after one time of flight T — N s /v g (= Tfl ight ). Since it grows 
linearly with system size N s , we plot G/N s . The data are 
obtained at p a — 0.75 and a = 0.25. 
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FIG. 11. The parked car correlator G(N a ,r) in the N 
phase; the garage is no longer macroscopically occupied, such 
that G does not scale with N a any more. The correlations 
light up, like a lighthouse beam, at every time-of-flight in- 
terval Tflig h t = T. The data are obtained at p D — 0.25 and 
a = 1. 
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FIG. 12. The parked car correlator G(N s ,t) at the 
C-N transition; it has a strong anticorrelation after one 
time-of-flight interval, T, even though for small r it decays 
linearly just like in the C phase (and with almost the same 
group velocity). The data are obtained at p = a = 0.25. 
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FIG. 13. Evolution of the parked car correlator G(N s ,t) 
through the phase boundary at fixed system size N s = 512 
along a line of constant group velocity: v g — 1 — 2p r with 
p r = 0.25; a — 1 in the C phase and at the C-N transition, 
while a = in the N phase. 
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FIG. 14. Schematics of soliton-type local perturbations 
diffusively spreading from the initial configurations (dashed 
lines) to one time-of-flight interval later (solid lines); (a) in 
the C phase, (b) at the C-N transition, and (c) in the N 
phase. 



0.06 , 

o.o5 ; 

0.04 
0.03 
0.02 
0.01 

o 



(a)32 
64 
128 
256 




(b)32 
64 
128 
256 



0.2 0.4 0.6 0.8 1 1.2 1.4 

t/N, z 



0.2 0.4 0.6 0.8 1 1.2 1.4 

T 1/Z /N B 

FIG. 15. The parked car correlator G(N s ,t) at the MC-N 
transition; the shape is very similar to that inside the MC 
phase, see Fig. 9. The data are obtained at p — 1/2 and 
a = 1. 
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